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Abstract 

The original Wigner function provides a way of representing in phase space 
the quantum states of systems with continuous degrees of freedom. Wigner 
functions have also been developed for discrete quantum systems, one popular 
version being defined on a 2N x 2N discrete phase space for a system with 
iV orthogonal states. Here we investigate an alternative class of discrete 
Wigner functions, in which the field of real numbers that labels the axes of 
continuous phase space is replaced by a finite field having N elements. There 
exists such a field if and only if N is a power of a prime; so our formulation 
can be applied directly only to systems for which the state-space dimension 
takes such a value. Though this condition may seem limiting, we note that 
any quantum computer based on qubits meets the condition and can thus be 
accommodated within our scheme. The geometry of our N x N phase space 
also leads naturally to a method of constructing a complete set of N + 1 
mutually unbiased bases for the state space. 



PACS numbers: 03.65.Ta, 02.10.-v 
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1 Introduction 



Given any pure or mixed state of a quantum system with continuous degrees 
of freedom, one can represent the state by its Wigner function [TJ |2] , a real 
function on phase space. The Wigner function acts in some respects like a 
probability distribution, but it differs from a probability distribution in that 
it can take negative values. The Wigner function has been widely used in 
semiclassical calculations, and it is also used to facilitate the visualization 
and tomographic reconstruction of quantum states. For a system with a 
single degree of freedom, one of the most interesting features of the Wigner 
function is this: if one integrates the function along any axis in the two- 
dimensional phase space — the axis can represent any linear combination of 
position and momentum — the result is the correct probability distribution 
for an observable associated with that axis [Hill]. 

Generalizations of the Wigner function have been proposed that apply 
to quantum systems with a finite number iV of orthogonal states, and the 
present paper continues this line of research. In 1974 Buot introduced a 
discrete Weyl transform which, when applied to a one-dimensional periodic 
lattice of N sites (with N odd), generates a Wigner function defined on a 
phase space consisting of an N x iV array of points 5J. Buot's work is related 
to earlier work by Schwinger jB] , who did not explicitly generalize the Wigner 
function but identified a complete basis of N 2 orthogonal unitary operators 
(elements of the generalized Pauli group — or discrete Weyl-Heisenberg group) 
that can be used to define an N x N phase space. A different approach was 
taken in 1980 by Hannay and Berry [7j: these authors directly adapted the 
definition of the continuous Wigner function to a periodic lattice and thereby 
arrived at a discrete Wigner function defined on a 2N x 2N phase space. 

Both of these basic approaches were later rediscovered and developed fur- 
ther by other researchers. Variations on the N x N scheme were proposed 
by Wootters jlj, Galetti and De Toledo Piza and Cohendet et al. 0, 
following initial investigations into the N = 2 case by Cohen and Scully 10 
and Feynman [TT]. 1 The N x N phase-space description has been applied to 

1 The Nx N approach has been problematic when N is even in that the method of Buot 
does not lead to a complete basis of Hermitian operators in that case (see Refs. P^l - 
and HH). In Ref. 0, the state-space dimension N is restricted to odd values; in Refs. 0] 
and |B] the difficulty is addressed by giving a special role to prime values of N. Schwinger 
likewise found it natural to regard each prime value of N as representing a single degree 
of freedom [§]. 
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quantum optics by Vaccaro and Pegg [15] and to quantum teleportation by 
Koniorczyk et al. [T^j- Discrete Wigner functions on the 2N x 2N model have 
been investigated by Leonhardt [T3| and used by Bianucci et al, Miquel et 
al., and Paz to analyze various quantum processes such as the Grover search 
algorithm JTj . All of these proposals have the feature that one can sum the 
Wigner function along different axes in the discrete phase space (including 
skew axes) to obtain correct probability distributions for observables associ- 
ated with those axes. Leonhardt in particular has emphasized the value of 
this feature for tomography, that is, for ascertaining the quantum state of 
a given ensemble by performing a series of measurements on subensembles. 
Other discrete Wigner functions have been considered which do not have this 
feature [TH] , but in the work we present here this tomographic property plays 
a central role. One can find further discussion of discrete Wigner functions 
and their history in, for example, Refs. [T2J EH EB] ■ 

In the continuous case, for a system with one degree of freedom, one can 
regard the Wigner function as being based on a certain quantum structure 
that one imposes on the classical phase space. The structure consists of as- 
signing to each straight line in phase space a particular quantum state. Let q 
and p be the phase-space coordinates, and suppose that the line in question 
is the solution to the linear equation aq + bp = c. Then the quantum state 
assigned to this line is the eigenstate of the operator aq + bp with eigenvalue 
c. Once this connection is made between lines in phase space and quantum 
states, one defines the Wigner function as W(q,p) = (l/2nh) Tr [pA(q,p)}, 
where p is the density matrix being represented and the operator A(q,p) is 
built in a symmetric way out of all the quantum states assigned to lines of 
phase space, the weight given to a particular state depending on the rela- 
tionship of its line to the point (q,p)- 

In this paper we wish to define a discrete Wigner function — actually a 
class of discrete Wigner functions — following as closely as possible the spirit 
of the construction just described. Because this construction is essentially 
geometrical, we want the geometry of our discrete phase space to be closely 
analogous to the geometry of an ordinary plane. For example, we need to 
have the concept of "parallel lines" in phase space, and we want two non- 
parallel lines always to intersect in exactly one point, just as in the Euclidean 
plane. Such considerations lead us to use, as the variables that label the axes 
of phase space, quantities that take values in a field in the algebraic sense. 
That is, for our axis variables q and p, we replace the usual real coordinates 
with coordinates taking values in Fjv, the finite field with N elements. (Our 
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phase space can therefore be pictured as an N x N lattice.) Now, there exists 
a field having exactly N elements if and only if N is a power of a prime 
[T§] . Thus our formulation is directly applicable only to quantum systems 
for which the dimension of the state space is such a number. It is always 
possible to extend it to other values of iV by taking Cartesian products of the 
basic phase spaces — the same strategy is used in Ref. jl], and indeed, exactly 
the same strategy is used in the continuous case when there is more than 
one degree of freedom — but in this paper we will restrict our attention to the 
basic phase spaces with field elements as coordinates. The use of arbitrary 
finite fields is what distinguishes our work from earlier approaches to discrete 
phase space. 

Though the restriction to powers of primes rules out many quantum sys- 
tems, there is one familiar case to which our formulation may be ideally 
suited, namely, a system of n qubits such as is commonly used to model a 
quantum computer. In that case the dimension of the state space is iV = 2 n , 
which is indeed a power of a prime. Thus our version of the discrete Wigner 
function provides an alternative to the 2N x 2N formulation that has been 
most frequently used in quantum information theoretic applications. Most 
likely each of these phase-space formulations will prove to have its own ad- 
vantages. 

As in the continuous case, we impose a quantum structure on the N x N 
phase space by assigning a quantum state to each line in phase space. We 
insist that this assignment satisfy a certain strong constraint, namely, that 
it transform in a particular way under translations. (The analogous quan- 
tum structure on the continuous phase space satisfies a similar constraint.) 
Any assignment of quantum states to lines that meets this condition we call 
a "quantum net," and we use it to define a discrete Wigner function. It 
turns out that the requirement of translational covariance does not pick out 
a unique assignment of quantum states to phase-space lines; that is, there is 
not a unique quantum net for a given N x N phase space. Moreover, we have 
not found a general principle that would select, in a natural way, one partic- 
ular quantum net for each N. So our approach does not lead immediately to 
a unique Wigner function for a quantum system with N orthogonal states. 
To some extent this non-uniqueness is mitigated by the fact that many dif- 
ferent quantum nets are closely related to each other. We define notions of 
"equivalence" and "similarity" for quantum nets and identify the similarity 
classes for N = 2, 3, and 4. A good portion of the paper is devoted to this 
classification of quantum nets, which amounts to a classification of possible 
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definitions of the Wigner function within this framework. 

One motivation for the present work comes from quantum tomography, 
which we mentioned above in connection with other discrete versions of the 
Wigner function as well as the continuous version. As we will see, our ap- 
proach leads naturally to a specific tomographic technique. Each complete 
set of parallel lines in the discrete phase space corresponds to a particu- 
lar measurement on the quantum system, or more precisely, to a particu- 
lar orthogonal basis for the state space. By experimentally determining the 
probabilities of the outcomes of this measurement, one can obtain some infor- 
mation about the Wigner function, namely, the sum of the Wigner function 
over each of those parallel lines. The sums over all the lines of phase space 
are sufficient to reconstruct the entire Wigner function and thus determine 
the state of the system. 

The particular orthogonal bases that are associated with sets of parallel 
lines turn out to be mutually unbiased, or mutually conjugate; that is, each 
vector in one of these bases is an equal-magnitude superposition of all the 
vectors in any of the other bases. Sets of mutually unbiased bases have been 
used before, not only for state determination [2H1 EE] but also for quantum 
cryptography and in other contexts [221 EH] , and a few methods have been 
found for generating such bases 1231 1201 I2ZI I2H1 E01 - As we 
will see, the discrete phase space developed in this paper leads to a rather 
elegant way of constructing mutually unbiased bases; it is essentially the 
same method as was discovered recently by Pittenger and Rubin j2H] and is 
closely related to the recent work of Durt [2H]) though those authors were 
not studying phase space or Wigner functions. The connection with mutually 
unbiased bases — valid for all prime power dimensions N — is one respect in 
which the Wigner function presented here is different from those proposed 
earlier. A consequence is that the tomographic scheme suggested by our 
phase space construction involves fewer distinct measurements than schemes 
derived from other discrete phase spaces (201 |2I]. This feature is the focus 
of Ref. [SI], which introduces for certain special cases some of the ideas that 
we present here in a more general setting. 

As further motivation, we note that the discrete Wigner function we de- 
velop here appears to bear an interesting relation to certain toy models of 
quantum mechanics proposed by Hardy [32J and Spekkens [33J to address 
foundational issues. For example, in both of these models a "toybit" has 
exactly four underlying ontic states, which could be taken to correspond to 
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the four points of our one-qubit phase space. 2 As has been suggested by 
Spekkens, the discrete Wigner function might therefore facilitate the com- 
parison between quantum mechanics and these toy theories [33J. 

Our discrete phase space is also related to some work on quantum er- 
ror correcting codes, which is similarly based on finite fields. (In Section 
4 we point out aspects of this relationship.) It is conceivable, then, that 
our Wigner function could be of particular value when representing certain 
encodings of quantum states. 

The remaining sections are organized as follows. Section 2 recalls the 
definition of the usual Wigner function and shows how it can be obtained 
from an assignment of quantum states to the lines of phase space. In Section 3 
we give the mathematical description of our discrete phase space and discuss 
its geometrical properties. Section 4 shows how to build a quantum net on 
this discrete phase space and shows that the bases associated with different 
sets of parallel lines must be mutually unbiased. The notion of a quantum 
net is then used in Section 5 to construct a discrete Wigner function. In 
Sections 6 and 7 we define our notions of equivalence and similarity between 
quantum nets and identify the similarity classes for small values of N. Finally 
in Section 8 we review our results and contrast the discrete and continuous 
cases. 

2 The Wigner function constructed from eigen- 

states of aq + bp 

Here we briefly derive the usual definition of the continuous Wigner function 
in a way that lends itself to generalization to the discrete case. The quantum 
system in question is a particle moving in one dimension, and the coordinates 
of phase space are the position q and momentum p. 

We begin by assigning a quantum state to each line in phase space. Con- 
sider the line specified by the equation aq + bp = c, where the real numbers 
a, b, and c are arbitrary except that a and b cannot both be zero. To this line 
we assign the unique eigenstate of the operator aq + bp that has eigenvalue 

2 For the case of a single qubit, our phase-space formulation is the same as in Refs. 0], 
[5] , and JT] , but it is already significantly different when one enlarges the system to a pair 
of qubits. 
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c. In the position representation we can write this operator as 



aq + bp = aq — ibh— , (1) 
dq 



and the relevant eigenstate \ip a bc) is given by 3 



The normalization of ip abc is chosen so that the integral \ip a b c ) (i> a bc\dc, 
with c 2 > Ci, is a projection operator. 

As we mentioned in the Introduction, given a density matrix p of the 
particle, the corresponding Wigner function will be of the form 

W(q,p) = ^TT[pA(q,p)}, (3) 

where A(q,p) is an operator that we will assign to the point (q,p). This 
operator is constructed, as we will see below, out of the states \ip a bc)- 

We want the Wigner function to have the property that its integral over 
the strip of phase space bounded by the lines aq + bp = C\ and aq + bp = c 2 
is the probability that the operator aq + bp will take a value between c\ and 
c 2 . This is one of the characteristic features of the Wigner function and is 
the property that makes it so useful for tomography. We can guarantee this 
property by insisting that the integral of (l/2irti)A(q,p) over the same strip 
of phase space is the projection operator onto the subspace corresponding to 
the eigenvalues of aq + bp lying between c\ and c 2 . That is, we insist that 

/ A(q,p)dqdp= f | ^ ahc ) (^abc\dc, for c 2 > c x . (4) 

J strip J ci 
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An equivalent expression of this condition, in terms of a single line in phase 
space rather than a strip, is the following: 

^ J 5(c-aq-bp)A(q,p)dqdp=\ip abc )(ip abc \. (5) 



3 For the special cases a = and b = we can take the eigenfunctions to be V'ofcc = 
(l/y/2irh\b\)e ic ^ hb and VaOc = (l/Va)S(q - c/a) respectively. 
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To find A(q,p) explicitly, we need to invert Eq. (J5J). But Eq. (JHJ) is an 
example of the well-studied Radon transform — the operator |Vw)(Vw|, re- 
garded as a function of a, b, and c, is the Radon transform of (l/2Trh)A(q,p) 
regarded as a function of q and p — and the inverse of this transform is well 
known [331 . Here we simply state the result: 



h\c\ f°° f°° ( 1 

A(l>P) = — / / ( ^ ( c - aq - bp) 2 ' ^ abc ^ abc ' da db ' ( 6 ) 
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where c is any nonzero real constant, and 1Z indicates the canonical regular- 
ization of the singular function that follows it. In the case of the function 
l/x 2 , this regularization is defined by 

74) /(,) * = r mutpzm ix . (7) 

Using the expression for ip a b c of Eq. (j2J), one can carry out the integration of 
Eq. © to get 

(q'\A(q,p)\q") = 5 (*±£ - q) e^^'X (8) 



The Wigner function then comes out to be 

1 1 r°° 

W(q lP ) = — Tr [ P A(q,p)] = — (q - x\p\q + x)e 2ixp/h dx. (9) 

Notice that according to the inverse Radon transform given in Eq. (|H)l. 
the operator A(q,p) is built out of all the operators \ip a bc) (ipabc\, the weight 
given to each operator tending to fall off as the associated line gets farther 
from the point (q,p). We will find that the analogous inversion in the discrete 
phase space is much simpler: the A operator associated with a given point 
is built entirely from the states assigned to the lines passing through that 
point. 

A particular property of the Wigner function that we want to generalize 
to the discrete case is translational covariance [2] . Here we state the property 
without proof. Let W(q, p) be the Wigner function corresponding to a density 
matrix p, and let p' be obtained from p by a displacement x in position and 
a boost y in momentum: 

pi — e i(y<i-xp)/hp e -i{y<i-xp)/f\ 
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Then the Wigner function W corresponding to p' is obtained from W via 
the transformation 

W(q,p) = W(q-x, P -y). (11) 

That is, when the density matrix is translated, the Wigner function follows 
along rigidly. 

Before moving on to the discrete phase space, let us mention an interesting 
property of the states l^w) that likewise has an analogue in the discrete case. 
Consider two infinite strips S and S' of phase space that are not parallel. The 
strip S is bounded by the lines aq + bp = c\ and aq + bp = c%, while 5" is 
bounded by a'q + b'p = c[ and a'q + b'p = c' 2 , and we assume that C2 > C\ 
and c' 2 > d x . Let P be the projection operator onto the subspace associated 
with S; i.e., 

P= I 2 \^a b c)(^a b c\dc. (12) 



Similarly, let P' be the projection onto the subspace associated with S' . 
Using Eq. (J2J) we can write down an explicit expression for P in the position 
representation (with a suitable modification if b = 0): 



(q\P\q') = - f 1 a sin 
n(q - q') 



(c 2 -ci)(g-g / ) 
2h\b\ 



e -(i/2hb)(q-q') [a(q+q')- (ci +c 2 )] / ]_ 3 ^ 



and P' can be written similarly. One can show by explicit integration that 
the quantity Tr(PP'), that is, f (q\P\q')(q'\P'\q) dqdq 1 , works out to be 

v ; 2nh \ab'-a'b\ V ; 

But the positive quantity {ci — ci)(c' 2 — c[)/\ab' — a'b\ is simply the area of 
the region where the two infinite strips overlap. Thus Tr(PP') is equal to 
this area expressed in units of Planck's constant. In the limit as the width 
of the strip S shrinks to zero, this result tells us that any eigenstate of the 
operator aq + bp yields a uniform distribution of the values of the operator 
a'q + b'p. 

As we will see, the analogue of this property in the case of discrete phase 
space is simpler. In place of strips we will consider individual lines of the 
discrete phase space. As we have said in the Introduction, each complete set 
of parallel lines will be associated with an orthogonal basis, and one finds 
that the magnitude of the inner product between any two vectors chosen 
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from different bases is always the same. This is the property called mutual 
unbiasedness. Before we can see how this comes about, and before we explore 
discrete generalizations of the Wigner function, we need to define our discrete 
phase space. 

3 Mathematical description of discrete phase 
space 

Our approach to generalizing the continuous phase space to the discrete case 
is quite simple. Like the continuous phase space for a system with one degree 
of freedom, our discrete phase space is a two-dimensional vector space, with 
points labeled by the ordered pair (q,p). But instead of being a vector space 
over the real numbers, it is a vector space over a finite field, and q and p are 
field elements. The number of elements in the finite field is the dimension N 
of the state space of the system we are describing. The physical interpretation 
of this discrete phase space will be left mostly to Section 4. In this section 
we focus on its mathematical properties. 

A field, in the algebraic sense, is an arithmetic system with addition 
and multiplication, such that the operations are commutative, associative, 
distributive, and invertible (except that there is no multiplicative inverse for 
the number zero) fH]- The real numbers are a familiar example of a field 
with an infinite number of elements. As we have said in the Introduction, 
there exists a field with exactly iV elements if and only if iV is a power of a 
prime, so our scheme applies directly only to quantum systems for which the 
state-space dimension is such a number. Moreover for any of these allowed 
values of N, there is essentially only one field having iV elements — any two 
representations are isomorphic — and we label this field F^v- If N is prime, 

consists of the numbers 0, 1, . . . , N — 1 with addition and multiplication 
mod N. If N = r n , with r prime and n an integer greater than 1, then the 
field Ftv is not modular in this sense but can be constructed from the prime 
field F r ; one says that F^r is an extension of F r . 

Let us illustrate this process of extension in the case of F 4 , which we 
will use frequently as an example. To generate F4, one begins by finding a 
polynomial of degree 2, with coefficients in F2, that cannot be factored in 
F 2 . (To generate F r n one would use a polynomial of degree n.) It happens 
that the only such polynomial is x 2 + x + 1 : there is no solution in F 2 to the 



10 



equation 

x 2 + x + l = 0. (15) 

The extension is created by introducing a new element uj that is defined to 
solve this equation, just as, in creating the complex numbers from the reals, 
one defines the imaginary element % to solve the equation x 2 + 1 = 0. Once 
uj is included, another element, uj + 1, is forced into existence, as it were, by 
the requirement that the field be closed under addition. One thus arrives at 
F 4 : 

F 4 = {0,l,w,u;+1}, (16) 

with arithmetic determined uniquely by the fact that uj satisfies Eq. (|15|). 
For example, we can square uj as follows: 

uj 2 = -uj - 1 = (-I)uj + (-1) = uj + 1, (17) 

where we have used the fact that —1 = 1 mod 2. Similarly, we have 

(uj + l) 2 = uj 2 + (1 + l)uj + 1 = uj 2 + 1 = (uj + 1) + 1 = uj (18) 

Following common practice we will frequently use the symbol uj to represent 
the field element uj + 1. The complete addition and multiplication tables for 
F4 are given here: 
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We now explore some of the geometric features of the phase space for a 
generic (but prime power) value of N. We picture the space as an iV x iV 
array of points (q,p), with q running along the horizontal axis and p along 
the vertical axis. For definiteness we place the origin, (q,p) = (0,0), at the 
lower left-hand corner. The phase space for N = 4 is shown in Fig. 1(a), and 
in Fig. 1(b) we show a possible physical interpretation of the axis variables 
if the space is being used to describe a pair of spin-1/2 particles. (The 
physical interpretation will be explained further in the following section.) 
We emphasize, however, that these pictures are not essential to our basic 
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Figure 1: The 4x4 phase space with axes labeled (a) by field elements, and 
(b) by spin states. 

construction. For example, we will often speak of a "vertical line," but this 
term is simply shorthand for a set of points of the form (q, y) where q is fixed 
and y can take any field value. 

More generally, a line in the N x N phase space is the set of points 
satisfying an equation of the form aq + bp = c, where a, b, and c are elements 
of Ftv with a and b not both zero. Two lines are parallel if they can be 
represented by equations having the same a and b but different values of 
c. Because the field operations are so well-behaved — especially since every 
nonzero element has a multiplicative inverse — the usual rules governing lines 
and parallel lines apply: (i) given any two distinct points, exactly one line 
contains both points; (ii) given a point a and a line A not containing a, there 
is exactly one line parallel to A that contains a; (iii) two lines that are not 
parallel intersect in exactly one point. Note that these propositions would 
not be true for general N if we were always using modular arithmetic, as has 
been pointed out in Ref. [14] . Consider, for example, the case N — 4. Under 
arithmetic mod 4 the points {(0, 0), (1, 2), (2, 0), (3, 2)} form a line, namely, 
the line that solves the equation p = 2q. But {(0,0), (1,0), (2,0), (3,0)} is 
also a line, and it shares two points with the first one. 

There are exactly N(N + 1) lines in our phase space, and these can be 
grouped into N + 1 sets of parallel lines. To see this, note that each nonzero 



12 



point (q,p) determines a line through the origin, namely, the line consisting of 
the points (sq, sp) where s takes all values in F^r. Let us refer to a line through 
the origin as a ray. Now, there are N 2 — 1 nonzero points, but each ray 
contains iV — 1 such points; so the number of rays is (N 2 — 1)/(N — 1) — N+l. 
Each of these rays then defines a set of N parallel lines. Let us call a complete 
set of parallel lines a "striation" of the phase space. 4 The five striations of 
the 4x4 phase space are shown in Fig. 2. One can observe there that the 
lines follow the three rules mentioned above. 



■ • • 



• • • • 



• • • • 



• • • • 



• • • • 



• • • ■ 



• • • • • • 

• • • • 



A • • • 
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Figure 2: The striations of the 4x4 phase space. 



Just as in the continuous case, one can speak of translations of the discrete 
phase space. A translation is the addition of a constant vector to each point 
of the space. For example, in the 4x4 phase space as pictured in Fig. 1, 

4 In Ref. 0] a similar set was called a "foliation," because in that case the elements 
of the set were sometimes higher-dimensional slices of a multi-dimensional space. Since 
the lines in our current construction are one-dimensional, they are more like "striae" than 
"folia." 
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translating by the vector (1, 0) has the effect of interchanging the first two 
columns and interchanging the last two columns. We will denote by T a 
the translation by the vector a. T a acts on points in phase space: T a f3 = 
a + (3. But we will also sometimes apply T a to an entire line, in which case 
it translates each point in the given line to yield another line (possibly the 
same as the original). 

Shortly we will need the concept of a basis for a field. A basis for the field 
F r n is an ordered set of field elements (ei, . . . , e n ) such that every element x 
in ¥ r n can be expressed in the form 

n 

x = s ^x i e i , (19) 

i 

where each Xi is in the prime field F r . There are typically many possible 
bases for a given field. In F 4 , for example, we could take as a basis, 

or (l,ou), or (a;, a;). Because in this paper we need to talk about bases for 
Hilbert spaces as well as bases for fields, we will often refer to the latter as 
"field bases." 

We will also need the concept of a dual basis, which in turn depends on 
the notion of the trace of a field element. The trace of a field element x is 
defined by 

tr x = x + x r + x r ' 2 H h x r "" 1 . (20) 

(We distinguish it from the trace of an operator by the lower case "tr".) 
Though this definition may seem quite opaque on first reading, the trace has 
remarkably simple properties, the most important for us being that (i) the 
trace is always an element of the prime field F r , (ii) tr (x + y) = tr x + try, 
and (iii) trax = atrx, where a is any element of F r . Now, given any basis 
E — (ei, . . . , e n ) for F r n, there is a unique basis E — (ei, . . . , e n ) such that 

tre^ = (21) 

where Sij is the Kronecker delta. This unique basis is called the dual basis of 
E |19] . We can immediately use the dual basis to obtain, for fixed basis E 
and field element x, the unique coefficients in the expansion (fT9*j) . Starting 
with that expansion, we multiply both sides by e~j and take the trace: 

n 

tr {xe~j) = Xjtr {t^j) = Xj. (22) 

i 

The expansion coefficients Xj will be used in the following section as we lay 
down a quantum structure on our discrete phase space. 
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4 Assigning a quantum state to each line in 
phase space 

We now need to supply our discrete phase space with a physical interpre- 
tation. We will do this by assigning to each line in phase space a specific 
pure quantum state as represented by a rank-1 projection operator. Let Q 
be the function that makes this assignment. That is, for each line A in phase 
space, Q(X) is the projection operator representing a pure quantum state. 
We will impose one condition on Q, translational covariance, to be defined 
shortly. A function Q satisfying translational covariance we will call a "quan- 
tum net." Later we will see how each possible choice of the function Q leads 
to a different definition of the discrete Wigner function. 

For N = r n where r is prime, our phase space applies most naturally to a 
system consisting of n objects (which we call "particles," though they could 
be anything), each having an r-dimensional state space. We assume that our 
system has this structure. 

We have seen in Eq. fllOj) the sense in which the continuous Wigner func- 
tion is translationally covariant. To define an analogous property in the dis- 
crete case, we need a discrete analogue of the unitary translation operators 

T( x ,y) = exp[%g - xp)/h\ (23) 

that appear in Eq. (jlUj) . That is, for each discrete phase-space translation 
T(x,y) with x and y in Fjv, we will define a corresponding unitary operator 
Tfay) that acts on the state space. In choosing these unitary operators, we 
are guided by the following considerations, (i) We want the multiplication 
of these unitary operators to mimic the composition of translations; that is, 
we insist that for any vectors a and /3 in phase space, 

T a T p « T a+f3 , (24) 

where the symbol ~ indicates equality up to a phase factor that might depend 
on a and (3. (The unitary operators of Eq. (}2*3*j) have exactly the same relation 
to the addition of continuous phase-space vectors.) (ii) There should be 
"basic" translations corresponding to unitary operators that act on just one 
particle. We make the connection between a translation vector (x, y) and 
individual particles by expanding x and y in field bases, allowing ourselves 
to use a different basis for each of the two dimensions of phase space. Thus 
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we write 

n 

x = ^ x e i€i (25) 
1=1 

and 

n 

y = Z>/*/i, ( 26 ) 

i=i 

where = (e 1; . . . , e n ) and F = (fx, ■ ■ ■ , f n ) are field bases, and we associate 
the coefficients x ei and y^ with the ith particle. (The symbols e and / 
are included in the subscripts to indicate which basis is being used in the 
expansion.) A translation that involves only coefficients having a particular 
value of i should be associated with a unitary operator that acts only on the 
ith particle. (Later we will discuss how much freedom we have in choosing 
the field bases E and F.) (iii) In the single-particle state space, we choose our 
unitary operators to be as analogous as possible to the continuous operators. 
Let (|0), . . . , \r — 1)) be a standard basis for the single-particle state space. 
Then in the space of the zth particle, the "unit horizontal translation," with 
x e i = 1 and yfi — 0, is associated with the unitary operator X defined by 

X\k) = \k + l) (27) 

with the addition being in F r , and the "unit vertical translation," with x e i = 
and — 1, is assigned the unitary operator Z defined by 

Z\k) = e 2wik/r \k). (28) 

The operators X and Z, which are generalized Pauli matrices introduced 
long ago by Weyl j^Hl, have been used by many authors in many contexts 
(often with non-prime values of r as in Ref. [SE1), including studies of discrete 
phase spaces IE] and mutually unbiased bases [2E1 12Z| • Except for phase 
factors, our general unitary translation operators are now fixed by Eqs. (|24|). 
f!27|) and (|28)1 . We write them as follows (and this equation fixes the choice 
of phase factors): 

T( x ,y) = X X ^Z V ^ ® • • - ® X Xm Z?t«. (29) 

We note that the operators Tr x , y ) play an important role in the theory 
of quantum error correction: they are normally taken as the basic error 
operators acting on an r n -dimensional state space (usually with r = 2). Of- 
ten the indices x e i and yfi labeling these error operators are treated simply 
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as elements of F r without the additional field structure that we have as- 
sumed. However, for some purposes it has been found useful to treat x and 
y as elements of the extension ¥ r n as we have done here. (See for example 

Refs. inn EU El mo 

In order for the Wigner function — defined later — to be translationally 
covariant, we insist that the function Q(X) be translationally covariant in 
the following sense: for each line A and each phase-space vector a, 

Q(T a X) = T a Q{X)Tl (30) 

That is, if we translate a line A in phase space by the vector a, the associated 
quantum state is transformed by T a . This is quite a strong requirement. To 
see why, consider the line consisting of the points (sx, sy) where x and y are 
fixed (and not both zero) and s ranges over the whole field Fjv- This line, 
and each of the lines of its striation, are all invariant under a translation by 
the vector (sx, sy) for any value of s. This means that in order to satisfy 
Eq. (|30p. the projections Q that we assign to these lines must commute with 
T(sx,sy) for each value of s. (If we were to represent the quantum states by 
state vectors rather than by projectors, the state vectors assigned to these 
lines would have to be eigenvectors of T/ S x,sy) f° r eac h value of s.) But this is 
impossible unless all of the operators Tr S x,sy) commute with each other. The 
basic operators X and Z obey the simple commutation relation 

ZX = rjXZ, (31) 

where 77 = e 2m ^ r ' . For each value of s, the unitary operator Tr sx sy \ can be 
written as 

T (s x,sy) = X^^Z^f 1 ® ■ ■ ■ <g> X {sx)en Z {sy) f n . (32) 

It follows from Eq. (|31|) that these iV operators commute with each other if 
and only if the following condition is met for all pairs of field elements s and 
*: 

n n 

J2( sx )ej(ty)fj = ^2{tx)ej{sy)fj, (33) 
i=i 3=1 

where the operations are those of F r ; that is, they are mod r. It turns out 
that this condition can be very simply expressed in terms of Fjy. We show 
in Appendix A that Eq. (}3*3*j) is satisfied for all values of x, y, s and t if and 
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only if the field bases E = (e 1; . . . , e n ) and F = (f%, . . . , f n ) are related by 
an equation of the form 

fi = we h i = l,...,n, (34) 

where w is any element of the field F^v- Thus, because we insist on transla- 
tional covariance, we are not free to choose the bases E and F arbitrarily. 
These bases enter into the definitions of the translation operators T a , and if 
the bases do not satisfy Eq. (|34p. there is no function Q(X) that is transla- 
tionally covariant with respect to these operators. 5 

Suppose now that E and F do satisfy Eq. so that the N operators 
T{sx,sy)i for fixed x and y and all values of s, commute with each other. These 
unitary operators are traceless and mutually orthogonal in the sense that 

Tr (v;,,, ;/: y;,,.,,,) =0 if s^t. (35) 

It follows that they define a unique basis of simultaneous eigenvectors (up to 
phase factors). 6 Thus as long as this condition on the field bases is satisfied, 
our requirement of translational covariance picks out a unique orthogonal 
Hilbert space basis to associate with each striation. (We will see shortly that 
translational covariance also requires Q to assign a different element of this 
basis to each line of the given striation.) Moreover, it follows from the work 
of Bandyopadhyay et al. [27] that these Hilbert space bases are all mutually 
unbiased. Specifically, Bandyopadhyay et al. show the following: if a set of 
N 2 — 1 traceless and mutually orthogonal N x N unitary matrices can be 
partitioned into N + 1 subsets of equal size, such that the N —1 operators in 
each subset are commuting, then the bases of eigenvectors defined by these 
subsets are mutually unbiased. Our operators T( x>y ) satisfy this hypothesis 
as long as the field bases E and F satisfy Eq. ([34)1 . 

5 As we have mentioned above, in some papers on quantum error correction the authors 
have indexed the error operators with elements of the field F r n . These authors have also 
insisted that the commutation relations among error operators be expressible in a simple 
way in terms of the field algebra [23 CM EH EH The condition Eq. does not seem to 
have appeared explicitly in these papers, but it may well be implicit. 

6 Here again there is a connection with the theory of quantum codes, in which one 
frequently considers sets of commuting error operators: a quantum stabilizer code is in 
fact a joint eigenspace of the operators of such a set |41l 1421 RTH I44| . However, in our 
case the commuting set is maximal, so that the subspace defined by a set of eigenvalues 
is spanned by a single vector. The vectors we define in this way are thus examples of 
stabilizer states. 
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Note that because there are N + l striations, the above argument — which 
again is closely related to Refs. [2E] and [2H| — shows that one can construct 
N + l mutually unbiased bases in N complex dimensions when N is a power 
of a prime. For a general value of N it is known [23 E] that the number of 
mutually unbiased bases cannot exceed N + l, and other papers have shown 
in other ways that this number is exactly iV+ 1 when N is a power of a prime 
[%m 1251 12H1 12Z| ■ Remarkably, the maximum number of such bases appears 
to be unknown for any value of N that is not a power of a prime (but see 
Refs. EDI 113 SHI EI which shed light on that problem). 

Let us find the Hilbert space bases that our construction assigns to the 
vertical and horizontal striations. The vertical lines are invariant under trans- 
lations by vectors of the form (0, s); so the Hilbert space basis associated with 
this striation consists of the simultaneous eigenvectors of the operators T( , s ). 
These operators take the form 

T ( o, s ) = Z s ^ ® ■ ■ ■ ® Z s f- (36) 

and are thus all diagonal in the standard basis; so their simultaneous eigen- 
vectors are simply the standard basis vectors 

\h) ® • •• <g> \k n ). (37) 

The horizontal lines are invariant under translations by the vectors (s, 0); so 
the Hilbert space basis associated with this striation consists of the simula- 
taneous eigenvectors of 

T( S) o) = X Sel <g> • • • <g) X Sen . (38) 
One finds that these vectors are 

\j n ), (39) 

where the single-particle states \j), notationally distinguished by the curved 
bracket, are given by 

ij)=^i> jfc ifc>- ( 4 °) 
v k=i 

(Again, i] = e 2m l r .) Note that for these two special striations — vertical and 
horizontal — the associated Hilbert space bases do not depend on the choice 
of field bases. This is typically not the case for other striations. 



19 



Let us see how this all works out for the case N = 4. First we arbitrarily 
choose (ei, e 2 ) = {00, 1) as the field basis for the horizontal translation variable 
x. One finds that the unique dual of this basis is (ei,e 2 ) = (l,o)). Thus in 
order to make translational covariance possible we should choose the field 
basis (/1, 72) for y to be either (l,u>) or some multiple thereof. We achieve a 
certain simplicity if we multiply by u to get (fx, f 2 ) = (00, 1). Then the 

basis for y is the same as the basis for x. Having made these choices, we can 
write down the unitary operator associated with any translation. Consider, 
for example, the following three vectors which are proportional to each other: 
(1, u>), (u, u), (o>, 1). In terms of our field bases, we can express these vectors 
as (0 ei + lea, l/i + 0/ 2 ), (lex + 0e 2 , lfr + l/ 2 ), (l ei + le 2 , 0/j + l/ 2 ). Thus 
according to Eq. ()29|) the unitary operators associated with translations by 
these vectors are, respectively, 



T M = Z ® X, T ( 



(0.1) 



(41) 



where in this case X and Z are the ordinary Pauli matrices, expressed in the 
standard basis as 



X 



1 

1 



1 
-1 



(42) 



One can verify that the three operators of Eq. (|41j) commute with each other. 
The unique basis of simultaneous eigenvectors is 



/ "I \ 



V » / 



/ 1 \ 

-1 

i 

\ 1 J 



( 1 \ 

1 

—i 
\ * J 



( 1 \ 



(43) 



This, then, is the basis that we associate with the striation containing the 
line {(0,0), (l,a;), (u, u>), (u, 1)}. 

In the same way we can figure out what Hilbert space basis is associated 
with each of the other striations. Fig. 3 shows the complete correspondence 
explicitly; each striation is labeled, in the left-hand column, by a point be- 
longing to the line in that striation that passes through the origin. The 
striations are listed in the same order as in Fig. 2. One can verify that these 
five orthonormal bases are mutually unbiased, as they must be. 

So far our construction only assigns a Hilbert-space basis to each stria- 
tion. Given our definition of T( x ,y) i n Eq. (|29jl . this assignment is completely 
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Figure 3: The five bases generated by the five striations. 



determined once we have chosen a field basis for each of the two dimensions 
of phase space. We now turn to the question of assigning a specific state Q(X) 
to each line A of phase space. How much freedom do we have in making this 
assignment? 

Consider a striation S. Let B = {\b s )} be the basis associated with this 
striation, with s e Fjv- We now consider a specific line Ago m S, namely, the 
ray that is included in S; that is, A50 is the line in S that passes through 
the origin. We are free to assign any of the states |6 S ) to \ so ; this choice is 
arbitrary. However, once we have made this choice, the vector assigned to 
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any other line of the striation is determined by Eq. (|30|): 

Q(T {x>y) \ S0 ) = T {x , y) Q(\ SQ )Tl y) , (44) 

since any line in the striation can be obtained by translating \$o- The func- 
tion Q(X) is thus entirely determined once we have assigned a quantum state 
to each of the rays of phase space. Moreover, it is clear from Eq. (}4*4*j) that the 
same quantum state cannot be assigned to two distinct lines of a striation: 
an operator Ti x y \ that translates A into A' cannot commute with Q(X), since 
it has a complete set of eigenvectors (not all degenerate) that are unbiased 
with respect to the basis associated with A. 

Let us summarize the choices we are allowed as we set up a quantum net 
for the N x N phase space. First, we choose a field basis E for the horizontal 
coordinate; any basis will do. Next, we choose a field basis F for the vertical 
coordinate, but here we are not so free. F must be a multiple of the unique 
basis dual to E; that is, fi = wti for some nonzero w G ¥ N . These choices 
determine the unitary translation operators according to Eq. (|29|h which in 
turn define a unique orthonormal basis to be associated with each striation. 
Now, for each striation S, we choose a particular vector \b) in that striation's 
basis and let Q(Xso) = \b)(b\, where Ag is the ray defining that striation. 
The state Q(X) assigned to any other line A is then determined uniquely by 
the condition Q(T (x ,y)X S o) = T( Xt y)Q(X S0 )T{ x>y y 

In the case N = 4, with the field bases E = F = 1) as before, we 
can define a quantum net by choosing, from each of the five bases shown 
in Fig. 3, one state vector to be associated with the corresponding ray in 
phase space. For example, we might choose, for each basis, the vector in the 
left-most column of that table. With this choice, the vertical line through 
the origin is associated with the state 1 00) [that is, k\ = k 2 = in Eq. (J3"7j) ]. 
and the other vertical lines, from left to right in Fig. 1(a), are associated 
with the states 1 01) , |10), and 1 11) respectively. If the system in question is 
a pair of spin-1/2 particles and if we interpret |0) as ||) and |1) as |f), the 
vertical lines can be labeled as shown in Fig. 1(b): ft, f|, |f, | L With this 
same choice, the horizontal lines are associated with the states | > >) and so 
on, as is also indicated in Fig. 1(b). 

In the next section we show how we can use a quantum net to define a 
Wigner function. 
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5 Defining a Wigner function 

A quantum net assigns a state Q(X) to each line A in phase space. The 
Wigner function W(q,p) of a quantum system should be such that when 
W(q,p) is summed over the line A, the result is the probability that the 
quantum system will be found in the state Q(X). That is, if p is the density 
matrix of the system, we insist that 



For a given quantum net Q, this condition completely determines the relation 
between p and W. 

We now use Eq. ()45|) to express W a explicitly in terms of p. We begin 
by observing that through any point a there are N + 1 lines, and that each 
point j3 7^ a lies on exactly one of these lines. These geometrical facts allow 
us to write 



where the first sum is over all lines that contain the point a. Using Eq. ([45 j) 
we can rewrite this as 




(45) 




(46) 




(47) 



where 




(48) 



Eq. pTjl is our explicit formula for W a . 

The operators A a have a number of special properties: 



1. A n is Hermitian. 



2. Tr A 



a 



1. 



3. TrA a Ap — N8 a p. 



4- E aeX A a = NQ(\). 
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These can all be proven directly from the definition. For our present purpose 
the most important is property (3), which we now prove explicitly. Starting 
with Eq. ([48)1 . we can write 

Tr A a Ap = J2Y1 ^WW)]- 2 E TrQ(A) +Tr/. (49) 

The last two terms have the value —2(N + 1) + N = —(N + 2). The value 
of the double sum over A and v depends on whether a and (3 are the same 
point. If they are, then of the (N + l) 2 terms in the sum, N + 1 of them have 
the value 1 because A = u, and the rest have the value l/N, because bases 
associated with different striations are mutually unbiased. Thus in this case 
we have 

TrA a Ap = (JV + 1) + N(N + 1)— - (N + 2) = N (a = (3). (50) 

If a ^ /3, then exactly one term in the double sum has the value 1, N terms 
have the value because A and v are parallel but different, and the rest have 
the value l/N because of the mutual unbiasedness. This gives us 

Tr A a Ap = 1 + N(N + 1)1 - (N + 2) = (a ^ (3), (51) 

which finishes the proof of property (3). 

Property (3) shows that the operators A a constitute a complete basis for 
the space of N x N matrices. In particular, we can write the density matrix 
as a linear combination 

p = J2b a A a , (52) 

a 

where the coefficients b a must be real since p and the A a 's are Hermitian. 
Multiplying both sides of Eq. by Ap, taking the trace, and using property 
(3) above, we find that b a is in fact equal to W a as expressed in Eq. p7j). 
We have thus found an explicit expression for the density matrix in terms of 
the Wigner function: 

a 

We now list a number of properties of the Wigner function and its rela- 
tionship to the density matrix. 

1. W n is real. 
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2. J2a&\ W a = Tr \pQ(\)]. This is the property (f43|) that we used to define 
the Wigner function. 

3. — 1- This follows immediately from property 2: break the sum 
over a into parallel lines, and the corresponding probabilities must sum 
to one. 

4. Let W be the Wigner function corresponding to a density matrix p and 
let W correspond to p', where p' = TppT^. Then W' a = W a _p. This is 
the translational covariance of the discrete Wigner function and is the 
analogue of Eq. ()10p. The proof is straightforward: 

W' a = ^Tr (p'A a ) = ^Tr {TppT^) 

= ^Tr { P TlA a T p ) = llr (pA a . p ) = W a ^ (54) 

Here we have used the fact that A a is a linear combination of the 
identity operator and the projections Q, which were constructed to be 
translationally covariant in accordance with Eq. QUO J). 

Of course the definition of the Wigner function depends on the quantum 
net Q; different choices of the quantum net will yield different definitions of 
the Wigner function. In order to show some examples of Wigner functions, 
for the remainder of this section we adopt the particular quantum net for 
N = 4 that we mentioned at the end of the preceding section. Recall that 
for this quantum net, we have taken the field bases to be E = F = (u, 1), 
and we have chosen to associate the rays of phase space with the states listed 
in the left-most column of Fig. 3. With these choices we can compute the 
operators A a and thereby find the Wigner function associated with any state 
p. In Fig. 4 we give the result for certain quantum states of a pair of spin-1/2 
particles, representing spin states as we did at the end of Section 4. 

One can check that the sum over any line is the correct probability of 
the state associated with that line. For example, in the case of the singlet 
state -75 ( I T 1) — lit)); if both particles are measured in the up-down basis, 
the only possible outcomes are fj. and J,f, corresponding to the two middle 
columns; similarly if both particles are measured in the right-left basis, the 
only possible outcomes are — > <— and < 

The property (|43jl of the discrete Wigner function is the one that makes it 
useful for tomography. Suppose that one has an ensemble of systems with an 
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Figure 4: Wigner functions for three states of a pair of qubits. 



iV-dimensional state space all prepared by the same process, so that each in- 
stance should be describable by the same (possibly mixed) quantum state. To 
find the values of the Wigner function, one can perform, on N + 1 subensem- 
bles, the orthogonal measurements associated with the N + 1 striations of 
phase space. From the probabilities of the outcomes one can reconstruct the 
Wigner function. In fact, from Eq. (J47|) one obtains the following equation 
for this reconstruction: 



1 

N 



\3a 



(55) 



where P(X) = Tr[pQ(X)] is the probability of the outcome associated with 
the line A. 

In this discussion we are assuming that N = r n where r is prime. A 
system with such a value of N can alternatively be described using the Wigner 
function of Ref. J|j, for which the phase space is the direct sum of n r x r 
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phase spaces. The tomography on this 2n-dimensional phase space requires 
(r + l) n different measurements, which is always greater than the number 
N + 1 = r n + 1 required by our present scheme. Indeed, for any value of 
N, N + 1 is the minimum number of orthogonal measurements needed to 
reconstruct a general quantum state, since a general density matrix contains 
iV 2 — 1 = (iV + 1)(N — 1) independent real parameters and each measurement 
provides only N — 1 independent probabilities. 7 

6 Classifying quantum nets 

According to our construction in Section 4, a quantum net is determined 
once we (i) specify the field basis for each of the two axes of phase space, 
and (ii) select, for each striation, a vector (from the basis associated with 
that striation) to be assigned to the line through the origin. For the purpose 
of this section, let us assume that the choice of field bases is fixed once and 
for all. We are still free to choose which vector to associate with each ray. 
How many possible quantum nets do these choices give us? The answer 
is N N+1 , since there are N + 1 striations, and for each one we can choose 
among N basis vectors. But these N N+1 quantum nets are not all greatly 
different from each other, and in some cases the definitions they generate of 
the Wigner function are closely related. In order to get a sense of the range 
of significantly different Wigner function definitions, we now begin to classify 
the possible quantum nets. For this purpose we define two relations between 
quantum nets: equivalence and similarity. 

Let us call two quantum nets equivalent if they differ only by a unitary 
transformation of the state space. That is, two quantum nets Q and Q' are 
equivalent if and only if there exists a unitary transformation U such that, 
for each line A, Q'(X) = UQ{X)lP. For example, Q' might be related to Q 
by a translation of the phase space, which by construction implies a unitary 
relation between Q and Q'. 

How many equivalence classes of quantum nets are there? To answer this 
question, note first that, regardless of what states a quantum net assigns to 
the vertical lines, because they are orthogonal — in fact they must be the basis 

7 Moving away from the simple tomographic model, there are many other schemes for the 
reconstruction of quantum states. In particular, one can use non-orthogonal measurements 
or adaptive measurements |48| , or one can perform arbitrarily many distinct measurements 

EH 
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states \ki) ® • • • ® \k n ) in some order — we can always find a unitary transfor- 
mation that will bring them to the same basis but in a standard order, the 
state |0) ® ■ ■ ■ <S> |0) being associated with the vertical ray. Moreover, we still 
have freedom, by a further unitary transformation, to change the phases of 
these states arbitrarily Thus the state assigned to the horizontal ray, a state 
that must already be one of the states |ji) <8> ■ ■ • <8> \j n ) [Eqs. and (|4L)|)]. 
can be brought, by changes in the phases of its components, to the particular 
state |0) (g> • • • ® |0). And this exhausts our unitary freedom. If two quantum 
nets, after having their vertical and horizontal states brought to a standard 
form in this way, are not now identical, then they must not have been equiv- 
alent to begin with, since there is no further unitary freedom. To find the 
number of equivalence classes, we simply have to consider the freedom that 
remains once the states associated with the vertical and horizontal lines are 
fixed. We still have N — 1 striations left, and for each one we still have N 
vectors that we can assign to the ray associated with that striation. Thus 
the number of equivalence classes is N 1 ^^ 1 . 

Note that the above argument also shows that if two quantum nets are 
equivalent, they must be related by a translation of the phase space. Start- 
ing with a given quantum net, one can generate N 2 equivalent quantum nets 
by translation, using the N 2 translation operators (including the identity). 
Thus each equivalence class must have at least N 2 elements. But since there 
are N N_1 equivalence classes and a total of N N+1 quantum nets, each equiv- 
alence class must have exactly N 2 elements, namely, the ones obtained by 
translation. 

In order to define the notion of similarity, we consider a different sort of 
transformation of the discrete phase space, namely, a linear transformation. 
That is, we imagine mapping each point a of phase space into a point a' = 
La, where L is linear over the field If we think of a as a column vector 
with components x and y, we can think of L as a 2 x 2 matrix with elements 
in the field: 




We call two quantum nets Q and Q' similar if and only if there exists a linear 
transformation L on the phase space, together with a unitary transformation 
Ul on the state space, such that for every line A, 

Q'(\) = UlQ(L\)U L . (57) 
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That is, Q' is unitarily equivalent not necessarily to Q itself but to Q act- 
ing on a linearly transformed phase space. A linear transformation can be 
regarded as a matter of changing the basis vectors of phase space, as a uni- 
tary transformation is a change of basis in the state space. In this sense two 
quantum nets are similar if they are related to each other by changes of basis 
in these two spaces. It turns out that Eq. (J57|) can hold only if L has unit 
determinant. For suppose that Eq. (J57j) holds for some L and Ul- Then from 
the fact that both Q and Q' must be translationally covariant [Eq. (|30|)]. it 
follows that for all phase space vectors a, 

U L T a Ul « T La . (58) 

In Appendix B we show that Eq. ()58|) can be satisfied for all a only if L has 
unit determinant. We show further that for every unit-determinant linear 
transformation L, there exists a unitary Ul such that Eq. ([58)1 holds. (See 
also Refs. EH] which address a different formulation of the same general 
problem.) 

This latter fact has an important consequence for classifying quantum 
nets. Given a quantum net Q(X), suppose that we construct another function 
Q'(A) = U j L Q{LX)U L , where L is a unit-determinant linear transformation 
and Ul is the unitary operator whose existence is guaranteed by Eq. ()58|). 
Then Q' is also a legitimate quantum net, translationally covariant with 
respect to the original translation operators T a . Thus any function obtained 
from a quantum net by linearly transforming all the lines of phase space, is 
itself a quantum net up to a unitary transformation. We will use this fact 
shortly in the classification of quantum nets. 

In the rest of this section we characterize the similarity classes of quantum 
nets for N = 2,3, and 4. For this purpose it is helpful to introduce a unitarily 
invariant function T of three phase-space points [I]: 

T aM = ^Tr(A a A^Ay) , (59) 

where A a is defined in Eq. (148(1 . Because T is not affected by a unitary trans- 
formation of the quantum net, it is constant over each equivalence class. 
Indeed, it follows from the orthogonality relation Tr A a Ap = N5 a p that the 
function V completely characterizes the quantum net up to a unitary trans- 
formation. Therefore, two quantum nets Q and Q' are similar if and only 
if the corresponding functions T and T' are related by a unit-determinant 
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linear transformation of the phase space, i.e., if r' a/37 = T LaL ^ L ^. Thus T can 
be used to distinguish different similarity classes. 

Similarity classes for N = 2 

For higher dimensions we will need to specify a field basis for each of the 
two phase-space dimensions, but in the case of a single qubit, there is no such 
choice, since the only field basis consists of the single number 1. The number 
of equivalence classes in this case is 2 2_1 = 2. To construct a representative of 
each one, we first fix Q(X) for the vertical and horizontal rays: to the vertical 
ray we assign the state |0) and to the horizontal ray we assign the state 
|0) = (l/v / 2)(|0) + |1)). As explained above, we have this freedom within an 
equivalence class. The only choice remaining then, which distinguishes the 
two equivalence classes, is the state to be assigned to the diagonal ray A^. 
This state must be one of the two eigenstates of XZ (that is, of o~ y ). Let 
us call these states \y + ) and \y~), defined by \y + ) = (l/y / 2)(|0) and 
= (l/v2)(|0) — z|l)). As it turns out, the two resulting quantum nets 
are similar to each other. To see the similarity, in Eq. ()57|) choose 

L=(J J) and U L =-L(\ _\). (60) 

One can verify that if we let Q(\d) = \y+){y+\, then after applying L and Ul 
as in Eq. (joT|) we obtain Q'(Xd) = \y-)(y- |, but the states assigned to the 
vertical and horizontal rays are unchanged. Thus there is only one similarity 
class for N = 2. 

Though we have not needed T to classify the similarity classes in this 
case, for comparison with other values of iV it will be helpful to see some 
of the values of this function. Here we give the values of r o 7 , where "0" 
indicates the origin and 7 is an arbitrary phase space point. The values of 
r o 7 are the same for both of the equivalence classes; so the following picture 
is valid for both. In this picture the value of roo 7 is written in the location 
defined by 7. (Recall that the lower left-hand corner is our origin, so the 
value written there is r oo-) 



(61) 



Here the factor 1/4 multiplies each term in the array. As must be the case, 
the two equivalence classes do differ in other values of T: when a, /3, and 
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7 are all different, T aj3l is complex, and the values for the two equivalence 
classes are related by complex conjugation. 

As we have seen, each quantum net yields a particular definition of the 
discrete Wigner function via Eq. ([47)1 . The fact that there is only one similar- 
ity class for N = 2 means, then, that there is essentially only one definition 
of the discrete Wigner function for N = 2 within the present framework. 
The allowed quantum nets differ from each other only by a rotation of the 
qubit (equivalence) and/or an antiunitary spin flip (similarity). Up to these 
modifications, the definition given in Eq. (|47|) for N = 2 agrees with the 
discrete Wigner function defined in Refs. fH\ |4"1 IHj. 

Similarity classes for N = 3 

For the three-element field there are two possible field bases: (1) and (2). 
Let us fix (1) as our field basis for each of the two phase-space dimensions. 
The number of equivalence classes of quantum nets for a single qutrit is 
3 3_1 = 9. Again we focus on a particular representative from each equivalence 
class by fixing the states assigned to the vertical and horizontal rays: to the 
vertical ray we assign the state |0), and to the horizontal ray we assign 
the state |0). The difference between equivalence classes then lies in the 
choices we make for the other two striations. The bases associated with 
these striations are 



(62) 



and 

1 |, (63) 
V 

where r] = e 27 ™/ 3 . We need to choose one vector from each of these bases to 
assign to the remaining two rays. We now use the values of r o 7 to help us 
identify the similarity classes. 

If we choose the first vector listed in each of Eqs. (Jo^j) and (JfiHj) , we get 
the following values of r o 7 (again, the position in the table indicates the 
value of 7): 









007 



1 


1 


1 


1 


1 


1 


1 


1 


1 



(64) 
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On the other hand, if we make any other choice, we find that the analogous 
table contains three zeroes lying along one of the lines of phase space. Here 
is an example: 
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1 



(65) 



Now, in the 3x3 phase space there are exactly eight lines that do not pass 
through the origin. Moreover, with a unit-determinant linear transformation 
acting on the phase space, we can move any of these eight lines into any other; 
thus, starting with the zeroes as in Eq. (J65j) . we can move them to any other 
such line. We saw earlier that if we modify a quantum net by applying a unit- 
determinant linear transformation to the phase space, the resulting function 
is, up to a unitary transformation, another quantum net. Therefore, as we 
use such transformations to move the zeroes among these eight lines, we are 
generating eight inequivalent quantum nets that by definition are in the same 
similarity class. We have thus accounted for all nine equivalence classes and 
have found that they lie in exactly two similarity classes: a class of eight as 
exemplified by Eq. (J63j) . and the special case shown in Eq. ()64|) which is in a 
similarity class by itself. 

Since there are two similarity classes for N = 3, there are also two quite 
different definitions of the discrete Wigner function. The simpler one, whose 
quantum net yields the T of Eq. (jp^j) . is the same as the one defined in 
Refs. jHIH]- The other one, with a T like that shown in Eq. ()65[). appears to 
be new. It necessarily has many of the features of the simpler definition — 
e.g., the sums of the Wigner function along the lines of any striation are 
the probabilities of the outcomes of a measurement associated with that 
striation — but it lacks some of the symmetry. It is not clear whether there 
is any physical context in which one would choose to use this less symmetric 
definition of the Wigner function. If there is, presumably it would be a 
context in which a particular quantum state, associated with the line along 
which r o 7 is zero, plays a favored role. 

Similarity classes for N = 4 

As always, we begin by fixing a pair of field bases for the two dimensions 
of phase space. For N = 4, let us adopt the bases we have used in our 
earlier example: we associate with each dimension the basis With 
the bases fixed, the number of equivalence classes of quantum nets in this 
case is 4 4_1 = 64. Referring to the list of bases in Fig. 3, we can generate 
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quantum nets from the 64 equivalence classes by choosing one state vector 
from each of the last three bases. To see how these 64 cases sort themselves 
into similarity classes, we again rely on r o 7 - Calculating r o 7 explicitly for 
various cases, one obtains many different arrays, among which the following 
four are representative: 
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Though the last two have comparable features, we note that it is not possible 
to change one of them into the other by a linear transformation of the phase 
space. The first two are clearly not related to the others or to each other by 
linear transformations since, for example, they have different values of r oo, 
which is invariant under linear transformations. 

The fact that these four arrays are not related by linear transformations 
shows that there are at least four similarity classes. In fact, by count- 
ing the number of different functions T a ^ that one can obtain by unit- 
determinant linear transformations (including the possibility of complex con- 
jugation, which does not show up in r o 7 ), one finds that the four examples 
illustrated above generate 64 distinct equivalence classes. We can conclude, 
then, that we have not left anything out and that there are exactly four 
similarity classes. 

Suppose that one has chosen one state vector from each of the five bases 
in Fig. 3, each vector being assigned to the appropriate ray of phase space. 
(Now we are not fixing a priori the vectors to be chosen from the first two 
bases.) It would be good to have a simple algorithm that would determine to 
which of the four similarity classes the resulting quantum net belongs. One 
could of course compute r o 7 for the given quantum net and compare the 
result with the arrays given in Eqs. (JBTIj) and ([67]). But in fact there exists a 
much simpler method, as we now explain. 

Let us label the four columns of Fig. 3 with elements of F 4 : from left 
to right, we label the columns with the values 0, 1, uj and uj. (This is 
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not an entirely arbitrary labeling. In writing down the bases in Fig. 3, we 
consistently used the same vertical translation operators to determine the 
order in each of the last four bases. The first basis cannot be obtained in 
this way and in that case we used the horizontal translation operators.) The 
column-labels can be used to specify which vector we have chosen from each 
basis: Let a be the label of the vector chosen from the first basis, h the label 
of the vector chosen from the second basis, and so on. For convenience we 
repeat in Fig. 5 the list of bases, with the new labeling scheme. Thus if h = uj, 
for example, the state vector chosen from the second basis (corresponding to 
the horizontal ray) is (1/2)(|00) + |01) - |10) - |11)). 

It turns out that there is a function D(a,b,c,d,e), taking values in F4, 
such that the value of D determines the similarity class of the quantum net 
defined by (a, b, c, d, e). In Appendix D we present a method for finding the 
function D. Here we simply state the result: 



D = uj(a + b + c) + uj ( a b c d e) 
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where we are using ordinary matrix multiplication to express the quadratic 
terms, all the operations being in F 4 . The correspondence between the value 
of D and the similarity class is as follows: The values D = and D = 1 
correspond, respectively, to the two similarity classes whose r o 7 arrays are 
shown in Eq. (|66j) : the values D = uj and D = uj likewise correspond to the 
similarity classes of Eq. (|57jl . 

To give an example, consider the specific quantum net we used earlier, 
obtained by choosing the first vector in each of the five bases. In this case 
a = b = c = d = e = and therefore D = 0; so the above correspondence 
predicts (correctly) that the quantum net obtained in this way is in the 
similarity class with Tooo — 25/16. 

If we adopt the convention of representing each equivalence class by the 
unique quantum net in that class that has a = b = 0, we obtain a simplified 
form of D: 

D = ujc + cd + ujce + ujde. (a = b = 0) (69) 

From this equation (or in other ways), one can easily determine the number 
of equivalence classes in each of the four similarity classes. One finds that 
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Figure 5: Labeling scheme 
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quantum nets for N — 4. 
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there are twenty values of the triple (c,d,e) for which D(c,d,e) = 0, so 
that there are twenty equivalence classes in the first similarity class shown in 
Eq. (|66J) . Similarly there are twenty equivalence classes in the other similarity 
class shown that equation and twelve in each of the two classes represented 
in Eq. (JBTj) . Thus the total number of equivalence classes comes out to be 
20 + 20 + 12 + 12 = 64, as it should. 

Similarity classes for larger N 

For larger values of N, it becomes more difficult to work out all the 
possibilities for the function r ai g 7 as we did above. We now outline another 
method for determining the number of similarity classes. 

We have seen that applying a unit-determinant linear transformation to 
a quantum net Q yields, up to unitary equivalence, another quantum net in 
the same similarity class. Thus we can regard the group of unit-determinant 
linear transformations as acting on the set of equivalence classes of quantum 
nets, and from this point of view the similarity classes are seen as the orbits of 
the group. According to a theorem in group theory, the number t of distinct 
orbits generated by a group G acting on a finite set is given by 

* = pE<^)> ( 7 °) 



geG 



where \G\ is the size of the group and <ft(g) is the number of elements in the 
set that are fixed by g G G. Since elements from the same conjugacy class 
fix the same number of elements, it is sufficient to calculate the number of 
quantum nets fixed (up to unitary equivalence) by one element from each 
conjugacy class and then multiply by the number of elements in that class. 
Using this method, one finds 8 that there are 11 similarity classes for N = 5. 

While we have not performed this calculation for higher values of N, we 
know that the identity always fixes all N N_1 equivalence classes of quantum 
nets, and one can show that the number of unit-determinant linear transfor- 
mations is exactly N 3 — iV; so the number of similarity classes must be at 
least 

N N-1 

— > N N ~ 4 (71) 

which grows very rapidly for large N. Therefore, within the current frame- 
work, if one is going to use a discrete Wigner function to describe, say, a 

8 Wc thank Robert Terchunian for pointing out that the number quoted in earlier drafts 
of this paper was incorrect and for computing the correct value. 
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large number of qubits, one has perhaps too many possible definitions of the 
Wigner function to choose from. Is there some further criterion that would 
naturally restrict the choice to, say, a single similarity class? 

When N is an odd prime, there always exists one similarity class with 
more than the required symmetry. We saw this above in the case N = 3, 
where for one of the similarity classes, r o 7 was independent of 7. In fact, 
whenever N is an odd prime, there exists a quantum net for which 

r^ 7 = lrr (aA ^ A7+7Aa) , (72) 

where (x, y) A (x', y') = xy' — yx' jlj. Indeed there is only one such quantum 
net up to unitary equivalence, as can be seen from the fact that every unit- 
determinant linear transformation leaves this particular T a ^ 7 unchanged. 9 So 
when N is an odd prime, there is one definition of the Wigner function (up to 
unitary equivalence) that stands out because of its high degree of symmetry. 

The sole similarily class for iV = 2 does not possess quite this degree of 
symmetry, but here one does not have the problem of too many possibilities. 

What if iV is a power of a prime? We have studied in detail only one such 
case, N = 4. In that case, of the 64 equivalence classes, it turns out that 
there are exactly two for which the matrix A a , defined in Eq. (|48jl. has the 
following special property: it is a tensor product of two single-qubit matrices. 
(For this condition it does not matter which point a we choose: if A a is a 
tensor product, then so is Ap, since the translation operator Tp_ a that relates 
them is itself a tensor product.) In the notation of Fig. 5, these two special 
equivalence classes are the ones for which, with a = b = 0, the triple (c, d, e) 
takes the values (0, 0, 0) and (a), u, 1). They are both in the same similarity 
class, since .D(0,0,0) = D(u,uj,1) = 0. Looking at the vectors in question, 
one sees that these two quantum nets are complex conjugates of each other. 

We can construct the A operators for these two special follows. 
Let A, ,, with x, y 6 {0, 1}, be the A operators derived from either of the 
quantum nets for N = 2. And let us express a point a in the 4x4 phase space 
as a = {x\uj + x 2 ,yiuj + 2/2), in which we are using our standard field bases 
for N = 4. Then one can show that the following two sets of tensor-product 

9 In arriving at Eq. 1)72(1 we have assumed that the field basis for the vertical axis 
(consisting of just one field element since N is prime) is the same as the basis for the 
horizontal axis. A different choice has the effect of multiplying the exponent by a constant 
factor. 
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operators correspond to quantum nets for N = 4: 




(£2,2/2) 



(73) 



and 




(£2,2/2)' 



(74) 



where the bar indicates complex conjugation. Moreover these two sets cor- 
respond to two distinct equivalence classes. 



In Ref. |3], Wigner functions for composite dimensions were constructed 
by taking tensor products of A operators for prime dimensions. We see now 
that at least for iV = 4, we can use this simple tensor-product construc- 
tion and at the same time produce a Wigner function with the tomographic 
properties defined by the lines of Ff. (That is, the tomography involves only 
N + 1 = 5 measurements rather than (r + l) 2 = 9 measurements.) It is 
interesting to ask whether something similar can be done for any power of a 
prime. This consideration might also be used to pick out one of the many 
possible definitions of the discrete Wigner function that our formulation al- 
lows for large N. But at present we do not know whether such tensor-product 
structures exist, within our current framework, for other powers of primes. 



So far in our classification of quantum nets we have been assuming fixed 
bases E and F in which to expand the phase-space coordinates q and p. We 
now ask how the range of possibilities expands when we consider all allowed 
choices of these bases. After the preceding discussion one might wonder why 
we would want to consider additional possibilities. Indeed for most practical 
purposes this is surely unnecessary, but for understanding the mathematical 
structure of our formulation, our classification scheme would be incomplete 
if we did not allow other field bases. 

Recall that we can choose any field basis E — (ei, . . . , e„) for the hori- 
zontal coordinate q. The basis for the coordinate p must then be of the form 
F = (fi, . . . , f n ) = (we±, . . . , we n ) for some field element w. What we want 
to know now is this: which of these choices lead to quantum nets that are 
not unitarily equivalent to the ones we have already discussed? 

The question is easily resolved. Suppose that we switch from one pair of 
field bases {E,F) to a different pair (E',F'). The effect of this switch is to 



7 Changing the field bases 
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change the translation operators from 



T ( g iP) = X qel Z p ^ ® • • • ® X qm Z p f 



(75) 



to 



(76) 



If there exists a unitary operator {7 such that for each point a 



(77) 



then given any quantum net Q'(X) based on the operators T^, we can define a 
corresponding quantum net Q(X) = WQ'(X)U whose translation properties 
are determined by the operators T a . Thus if Eq. (|77|) is satisfied for some 
U, the change of field bases has not produced any new quantum nets, up to 
unitary equivalence. Now, we can identify two elementary kinds of change 
in the field bases that are allowed by the condition fj = wef (i) change ej 
arbitrarily into e'-, and simultaneously change fj into /• = we' j (with the 
same w as before); (ii) leave unchanged and change fj into /j = w'tj. Any 
allowed change of the field bases can be regarded as a combination of these 
two. Appendix C shows that under a change of the first kind, there exists a 
unitary operator U such that Eq. (J77)l is satisfied. Thus these changes do not 
produce any new equivalence classes of quantum nets. On the other hand, if 
we make a change of the second kind, we can write the resulting X" as 



Except in the trivial case where we have made no change at all, the deter- 
minant of this matrix is not unity, and therefore, as shown in Appendix B, 
there exists no unitary U such that Eq. (|77|) is satisfied. Thus this second 
kind of change of basis does produce new quantum nets. By performing such 
basis changes, we can multiply by iV — 1 the number of equivalence classes 
of quantum nets, since there are N — 1 choices for the non-zero field element 
w. 

In Fig. 6 we summarize in tabular form our classification of quantum 
nets for iV = 2, 3, and 4. Each box in the figure represents a similarity class, 
and the integer appearing inside the box indicates the number of distinct 
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equivalence classes within the given similarity class. The similarity classes 
are arranged in columns corresponding to different values of the field element 
w that expresses the relation between the bases E and F. Thus, for example, 
there are altogether 192 distinct equivalence classes for N = 4. In general 
the number of equivalence classes, now that we are allowing alternative field 
bases, is (N - l)N N -\ 

w = 1 

N = 2: 
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w — 1 w = 2 
N = 3: fT] [T] 

m m 
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Figure 6: Classification of quantum nets for iV = 2, 3, and 4. 



8 Discussion 

The main new contribution of this paper has been to use the general concept 
of a finite field to construct discrete phase spaces, and to study generaliza- 
tions of the Wigner function defined on such spaces. In this formulation, 
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there is not a unique definition of the discrete Wigner function for a given 
system; rather, the definition depends on the particular quantum structure 
that one lays down on the discrete phase space. This quantum structure, 
which we have called a quantum net, assigns a pure quantum state to each 
line in phase space. The assignment is severely constrained by the condition 
of translational covariance, which is analogous to a similar property of the 
continuous Wigner function. In particular, the quantum states assigned to 
parallel lines are forced by this condition to be orthogonal, and the orthogo- 
nal bases assigned to distinct sets of parallel lines are forced to be mutually 
unbiased. Because of this, our construction provides a method (closely re- 
lated to the methods of Refs. [2Hj and j2H]) of generating complete sets of 
mutually unbiased bases. 

It is interesting to contrast the discrete Wigner functions presented in 
this paper with the usual continuous Wigner function. In addition to trans- 
lational covariance, the usual Wigner function has another remarkable prop- 
erty which can be called covariance with respect to unit-determinant linear 
transformations |3H I5U] . Let p be any density matrix for a system with one 
continuous degree of freedom, and let W p (a) be its Wigner function, where 
a — (QiP) is a phase-space point. Now consider any unit-determinant lin- 
ear transformation L acting on phase space. It is a fact that for any such 
L, there exists a unitary operator Ul such that W p (La) = W p '(a), where 
p' = UlpU^. In other words, rotating the phase space, or stretching it in 
one direction while squeezing it in another by the same factor, is equivalent 
to performing a unitary transformation on the quantum state. That is, this 
sort of transformation of the Wigner function can in principle be carried out 
physically. The analogous property typically does not hold for our discrete 
Wigner functions. We can see this even in the case N = 2. In that case the 
linear transformation 



interchanges horizontal lines with vertical lines while leaving the diagonal 
lines unchanged. For any of our quantum nets, this corresponds to an inter- 
change between eigenstates of X and eigenstates of Z, while the eigenstates 
of XZ (or of Oy) remain unchanged. No unitary operator can effect such a 
transformation; so this L cannot be realized physically. 

Note that in our formulation one does find a weaker version of this prop- 
erty. Every unit-determinant linear transformation, while not necessarily 
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corresponding to a unitary transformation of the quantum state, does corre- 
spond to a unitary transformation, up to a phase factor, of the translation 
operators, as is shown in Appendix B. Moreover, there are certain special 
quantum nets for which the associated Wigner function does in fact have the 
stronger property These are the quantum nets discussed in Section El with 
T given by Eq. (J72J). But such special quantum nets appear to exist only 
for odd prime values of N. If one wants to generalize the Wigner function 
to other finite fields, including even the case of a single qubit, evidently one 
must do without some of the symmetry of the continuous Wigner function. 

There is another interesting difference between the continuous case and 
the discrete case. It is central to our construction that every line of dis- 
crete phase space corresponds to a quantum state, as is also true for the 
continuous phase space. However, in the continuous case, there is a specific 
correspondence between lines and quantum states that arises naturally: the 
quantum state assigned to the line defined by aq + bp = c is precisely the 
eigenstate of aq + bp with eigenvalue c. This correspondence is possible in 
part because the parameters a, b, and c used in the equation for the line also 
make sense as coefficients in the algebra of operators. In the discrete case, 
on the other hand, the parameters a, b, and c are elements of a finite field 
and cannot be combined in the same way with operators on a complex vector 
space. This is why, in the discrete case, there is not a unique quantum net 
for a given phase space. The requirement of translational covariance forces a 
certain correspondence between striations and bases, but not between lines 
and state vectors. 

In this connection, it is interesting to ask what new possibilities would 
open up in the continuous case if one were to approach the construction 
of distribution functions on continuous phase space along the lines we have 
followed in this paper. That is, rather than adopting a priori a particular 
correspondence between lines and quantum states, suppose that we were 
to allow, for each striation, a separate translation of the quantum states 
assigned to that striation. Most of the "generalized Wigner functions" that 
would thereby be allowed would no doubt be quite ugly, but one can imagine 
certain special quantum nets with useful properties. 

At one level what we have been exploring in this paper is the general 
concept of phase space. This concept is certainly central to the physics of 
systems with continuous coordinates. Just as certainly, it has been less cen- 
tral to the physics of discrete systems. However, as we have seen, even in 
the discrete case the notion of phase space, with axis variables taking values 
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in a field, meshes nicely with the complex-vector-space structure of quantum 
mechanics. The sets of parallel lines in phase space correspond perfectly with 
a complete set of mutually unbiased bases for the state space, and transla- 
tions in phase space correspond to physically realizable transformations of 
quantum states. Indeed, if one were starting with the complex vector space 
and the concept of mutually unbiased bases, and were trying to find a com- 
pact way of expressing quantum states in terms of such bases, one might be 
led naturally to phase space as the most economical framework in which to 
achieve this expression. 

At present we have no particular evidence that discrete phase space holds 
as distinguished a place with respect to the laws of physics as continuous 
phase space does. On the other hand, as a practical matter discrete phase 
space descriptions have been found useful in a variety of problems in physics 
(see for example Refs. [S], [7j, [TE], [H] and pTT]). and we hope that our 
phase space based on finite fields will find similar applications, especially 
in analyzing systems of qubits. Indeed, our formulation (as presented in a 
preprint) has already been applied by Galvao to a question regarding pure- 
state quantum computation [52] • Galvao makes explicit use of the full range 
of definitions of the Wigner function that our scheme allows. For other 
applications, it is likely that further research will have to be done to identify, 
out of the set of possible Wigner functions, a much smaller number in which 
the processes of interest are most simply represented. 
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Appendix A: Commuting translation opera- 
tors and the choice of field bases 

Recall the necessary and sufficient condition (J3*3*j) for the commutation of 
translation operators corresponding to parallel translations: 




(80) 



j 



3 
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Here we show that this condition is true for all s and t, and for all (x, y) 7^ 
(0,0), if and only if the field bases satisfy fj = wej for some nonzero field 
element w. 

We begin by assuming that fj = wej and proving that Eq. (JHUj) follows. 
From fj = wej it follows that tr(eifjW~ l ) = Sij. Thus we can write 

J2ji sx )ej(ty)fj = J2ii( sx )ei(ty)fjtr( y e i f j w- 1 ) 



(81) 



= tx(sxtyw~ 1 ) = tiitxsyw^ 1 ) 

which proves Eq. ([80]). 

Now we go the other direction. Assume Eq. ()80|) and note that for any 
z G Fjy, z e j = ti{e jz). Thus 



^2tr(ejSx)tr(fjty) = tr(e i tx)tr(/ j sy). 



(82) 



Using the linearity of the trace, we can rewrite this equation as 



tr 



sx 



tr 



Cjtx 



For this to be true for all x, we must have 



J2 tic (fjty)^ 



^ti(fjsy)e 3 



It follows that the quotient 



£j ^(fjty)ej 



is independent of t, though it might depend on y. That is, 

^2^(fjty)ej = A y t 



(83) 



(84) 



(85) 



(86) 



44 



for some A y e But that this is true for all t implies that 

fjU = Aytj for all y and j, (87) 

which in turn implies that A y is a constant times y and that fj is a constant 
times ej. Finally, the latter condition is equivalent to 

fj = wej, (88) 
which is what we wanted to prove. 

Appendix B: Changes in the translation oper- 
ators due to linear transformations 

1. A linear transformation that preserves the translation operators up to a 
unitary transformation and phase factors must have unit determinant. 
Recall the definition of the translation operators: 

r (w ,) = X qi Z Pl ® • ■ ■ <g> X q "Z Pn , (89) 

where we have suppressed the notation indicating the field bases in which q 
and p are expanded, since we are not going to be changing the bases in this 
section. 

Let L be a linear transformation of the phase space, and suppose that 
there exists a unitary operator Ul such that for every point a, 

U L T a Ul = eW^TLa, (90) 

where is any real function of L and a. We show now that this can be the 
case only if L has unit determinant. 

Consider the operator T a TpT%rl. Using the fact that ZX = rjXZ, where 
i] = exp(27rz/r), and the fact that Z r = X r = I, one finds that this operator 
simplifies to 

T a T p TlTl = 7/-™)/, (91) 

where a = (q,p), (3 = (x,y), and the dot product x ■ p stands for Yli x iPi- 
Now, if Eq. (}9Ti|) is true, it follows that 

U L T a T p TMul = T La T Lfj Tl a Tl , (92) 
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and since T a TpT£Tp is proportional to the identity operator, we can say 

T a T p TlTl = T La T Lp Tl a Tl p . (93) 

From this and Eq. (J91|) it follows that L must preserve the quantity (x-p—q-y), 
regarded as an element of F r . 

Let us now invoke the field bases E and F, which must satisfy the con- 
dition fj = wtj for some field element w. Since tr^ej) = 5^, we can write 

x ■ p — q ■ y = tr[w (xp — qy)}. (94) 

Thus the latter quantity must be conserved by L. Now, when (q,p) and 
(x, y) are both transformed by L, the effect on xp — qy is multiplication by 
the factor detL. So tr[u; -1 (:cp — qy)] must equal ti[w~ l (xp — qy) detL] for 
every (q,p) and (x,y). This is the same as saying that tr(6 detL) = tr(6) 
for every field element b, which is true only if detL = 1. Thus any linear 
transformation for which Eq. (jDUj) is valid must have unit determinant. 

2. Every unit- determinant L preserves the translation operators up to a uni- 
tary transformation and phase factors: the case N = r n with r an odd prime. 
Here we show that for every unit-determinant linear transformation L on 
phase space, there exists a unitary transformation Ul on state space such 
that for every phase-space vector a, 

U L T a Ul w T La , (95) 

restricting our attention for now to the case where N is a power of an odd 
prime. In this case we can, as we see below, specify the phase factor that is 
implicit in Eq. (|93j): 

U L T a U{ = rjMMtf-riT^ ( 96 ) 

Here 77 = e 2nl ^ r , (q,p) = a, and (q',p') = La. Also, in Eq. (|96|) and in what 
follows, the exponent is first computed as an element of the field F r and is 
then interpreted as an integer in the set {0,...,r — 1}. (For example, if 
r = 3, the expression i] 1 ^ 2 is interpreted as 7] 2 , since in F3, 1/2 has the value 
2.) To prove that such a Ul exists, we define the following linear mapping 
M on the space of N x iV matrices: 

M(T a ) = r}W 2 M- p '- q -rtT La . (97) 
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This equation defines M on all the translation operators and thus by linearity 
on all operators. Our aim is to show that the M defined by Eq. (|97|) is of the 
form M(B) = UBU^ for some unitary operator U. We do this by showing 
first that M preserves multiplication; that is, for any phase-space vectors a 
and /?, 

M(T a T (3 )=M(T a )M(T f3 ). (98) 
This we do by direct calculation, starting with 

T a T p = (X qi Z Pl ®---®X qn Z Pn )(X Xl Z yi ®---®X Xn Z Vn ) =ri ix - p) T a+ p. (99) 

Here (x, y) = j3, and we have used the fact that 

ZX = rjXZ. (100) 

Thus 

M(T a Tp) = ^y^+^p'+v')-^)<P+y)^^p)T L{a+ py (101) 
On the other hand, M(T a )M(Tp) is given by 

= ?? (i/ 2 )( g 'y- g . P ) r/ (i/2)(.'v-,. ?/ ) r/ ( a: 'y) Ti(Q+/3) . (102 ) 

Comparing the exponents in Eq. (|1U1J1 and (|102|) . one finds that they are 
equal (as elements of F r .) as long as q' ■ y' — x' ■ p' is equal to q ■ y — x ■ p. 
But this condition is guaranteed by the fact that L has unit determinant [see 
Eq. (|94|)]. So M does indeed preserve multiplication. 

A linear transformation on the set of all N x N matrices that preserves 
multiplication must be a conjugation; that is, there must exist a matrix S 
such that for any N x N matrix B, 

M(B) = SBS-\ (103) 

But M has another special property that we now prove, namely, that for any 
N x N matrix B, 

M(B ] ) = [M(B)]l (104) 

Let us show that Eq. ()104j) is satisfied when B is any of the translation 
operators T a ; it will then follow that the equation is true for any B. First, 
we have 

Tt = Z~ Pl X- qi <g> • ■ • ® Z- pn X~ q " = rj {q - p) T_ a . (105) 
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Thus 



M(Tl) = V (q ^M{T_ a ) = v (i-p) v (mW-P'-g-P) T _ La = ^ l ' 2 ^ p+q '^T_ La . 

(106) 

And 

= r/^^'+^T^. (107) 
From Eqs. f)106j) and ()107|) we see that 

M(Tl) = [M(T a )]\ (108) 

from which Eq. ()104J) follows. But Eq. (|104|) cannot be true for all B unless 
the matrix S in Eq. (|103|) is unitary. This proves the desired result when N 
is a power of an odd prime. 

3. Every unit- determinant L preserves the translation operators up to a uni- 
tary transformation and phase factors: the case N = 2 n . 
The above proof does not work when N is a power of 2, because the division 
by 2 that appears in many of the exponents cannot be done in F2«. For 
this case we explicitly construct the desired unitary transformation, which 
we imagine acting on a system of n qubits. Since the case N = 2 n is the one 
most likely to be relevant for quantum computation, our explicit construction 
may also have a practical value. 

One can show that the group of unit-determinant linear transformations 
on F| n can be generated by the following three elements of the group: 

*-(!:) *-(!!) (109) 

where z is any primitive element of the multiplicative group of F 2 «; that 
is, any nonzero element of can be written as a power of z. (Such a z 
exists for any finite field [TUj.) Our plan is first to choose a specific pair of 
field bases E and F for the two phase-space coordinates, and then to find a 
unitary Ui for each L iy such that 

U t T a V\^T Ua (110) 
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when T a is denned in the chosen bases. We will argue separately that this 
result survives changes in the field bases. The specific bases we choose for now 
are the following: for the horizontal coordinate we use E = (e±, . . . ,e n ) = 
(1, z, z 2 , . . . , £ n_1 ), which is indeed a basis as long as z is a primitive element, 
and for the vertical coordinate we use the dual of E. 

Consider first the transformation L\. Acting on a generic phase-space 
point a = (q,p), it yields 



1 \ / q \ = ( q 
1 1 A P J V 1 + P 



[111) 



so that U\ must effect the transformation (up to a phase factor) 

X qi Z Pl ®---®X qn Z Pn -> X qi Z qi+Pl <8) ... (8) X qn Z qn+Pn . (112) 

Because the components qi and pi are not shuffled among the various qubits 
in this case, it is not hard to find a suitable Ui, the following is one of a 
number of operators that would suffice: 

*-(S!)--(i!)- ("3) 

If we think of the qubits as spin- 1/2 particles, this operator rotates each 
qubit by 90° around the z axis. Similarly, one finds that Ui can be taken as 
a tensor product, 

which rotates each spin by 90° around the x axis. Note that these definitions 
of Ui and U2 would have the desired effect regardless of the field bases we 
were using to expand q and p. 

We now consider the transformation L 3 . Our operator U 3 , which does de- 
pend on the chosen field bases, is constructed from two basic gates: CNOTij 
acts on qubits i and j, taking \k iy kj) to \k iy kj + fcj) with each index k G F 2 ; 
and SWAPij interchanges qubits i and j. In terms of these gates, U 3 is 

U 3 = (f[ CNOT^ SWAP ln SWAP 1{n . 1) ■ ■ ■ SWAP 12 . (115) 
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Here the a^s are the coefficients in the expansion of z n in the basis (1, z, z 2 , . . . , z n 1 ): 

n 

z n = J2a ] z 1 - 1 . (116) 

We note for future reference that a x must be equal to 1; if it were not, we 
could divide both sides of Eq. ()116|) by z and conclude that 1, z, z 2 , . . . , z n ~ x 
are not linearly independent, contradicting the fact that these elements form 
a basis. 

To see that U3 has the desired effect, it is sufficient to check its action on 
the basic translation operators 7(1,0), 7( z ,o), ■ ■ TLn-i ) and 7'( j 1 ), Tqj^, 
. . . , T(o,/„), where the /j's constitute the dual basis. In the first of these sets, 
consider for example 7( 1)0 ). Applying L 3 to the point (1,0) gives (z,0), so 
that we want U3 to effect the transformation 

X®/®---®/ -> I®X®/®---<g)I. (117) 

The Z/3 of Eq. ()115|) does accomplish this shift through the SWAP opera- 
tions; the CNOTs have no effect since by the time they act, the operator in 
the first position is the identity. The same sort of shifting operation works 
also for the other basic horizontal translation operators, with the exception 
of T( z n-i )0 ). In this last case, since L 3 takes (z n_1 ,0) to (z n ,0), we want U3 
to have the following effect: 

7®---®/®X -> X ai ®---®X an . (118) 

One can verify that this is indeed the effect of the U3 defined in Eq. (I115|) : 
now the SWAPs shift the operator X from the nth position to the first 
position, and the CNOTs change the operator J to X in every position j 
for which aj = 1. Here we also need the fact, mentioned above, that ai = 1, 
since the CNOTs will not affect the operator X in the first position. 

Before we consider the vertical translation operators Tr jA, it is helpful 
to introduce a matrix representation of multiplication by z. Let the matrix 
z, with components in F 2 , be defined by 

zej = Zjitj. (119) 
3 

One can show that, for any field basis E, the effect of z on the dual basis 
can be expressed as 

zej = Zjjej. (120) 
3 
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That is, one uses the transpose of the original matrix. For our particular 
basis, we have 

zet = e i+ i, i = l,...,n-l (121) 

and 

n 

ze n = ^a j e j (122) 
i=i 

It follows, then, that the effect of z on the dual basis is given by 

zfi = fn (123) 

and 

zfi = f-i + difn, i = 2,...,n (124) 

Now, because L 3 multiplies the vertical coordinate by z~ l , we want U3 to 
take T(oj.) to TV^-ijy (up to a phase factor) for each value of i. This means 
that U\ should take T^j^ to T( z j.). Eqs. f|123|) and ()124|) tell us how to 
write zfi as a sum of basis elements. And U\ is given by Eq. (jl!5J) but with 
the operators in the reverse order. By comparing the effect of U\ with the 
effect of multiplication by z, one can check that Jj\ does indeed transform 
T(0, fi) as desired. Since every translation operator can be written as a 
product of the basic horizontal and vertical translation operators, it follows 
that Eq. ()110|) holds for every T a . 

So far we have restricted our attention to translation operators defined 
in terms of a particular pair of bases. The following Appendix shows that 
Eq. (jimp can be extended to any pair of bases E and F, as long as f = wei 
for some field element w. 



Appendix C: The effects of changes in the field 
bases 

1. Changing from e-i and wci to gi and wcji 

The translation operators depend on the choice of two bases for the field, 
one for each coordinate. Here we ask how the translation operators change 
when we make a change of the following form in these two bases. Let the 
initial bases be E = (ei, . . . , e n ) and F — (fi, ■ ■ ■ , f n ) — ( w &i, ■ ■ ■ , wen). The 
translation operators in these bases are 

T a = X^ZPf 1 ® • • • ® X qen Z p f", (125) 



51 



where (q,p) = a. We now change the bases to G — (gi, . . . ,g n ) and H = 
(hi, . . . , h n ) = (wgi, . . . , wg n ), where G is an arbitrary basis and w is the 
same field element as in the definition of F. The translation operators arising 
from these bases are 

T' a = X qgl Z Phl <g> ■ ■ ■ <g) X qan Z Phn . (126) 

We show here that there exists a unitary operator U such that for every point 
a, 

T' a = UT a Ul (127) 

Our method is the same as in part 2 of Appendix B. We define a linear 
map M such that for every a, 

M(T a ) = T' a , (128) 

the action of M on other matrices being determined by linearity. We show 
that M preserves matrix multiplication and the adjoint operation and must 
therefore be conjugation by a unitary operator. 

First we look at the relation between the components of the same field 
element in two different bases. One can show that 

'/.</< zZ' : ^r ( 129 ) 

3 

where 7^ = tr(^ej). Let z/^ be the matrix that similarly expresses the 
relation between F and H: 

Phi = ^2vijPfj- (130) 

3 

We can see that v and 7 are closely related: 

u i:j = tr(hifj) = tT[(w~ 1 g i )(we j )] = tr(^ej) = (l~ r )ji- (131) 

That is, v is the transpose of the inverse of 7. 

We now show that for any a and f3, M(T a Tp) = M(T a )M(Tp). As we 
saw in Appendix B, 

T a T p = rf e ' P! T a+ p, (132) 

where (x, y) = 13. Thus 

M(T a T p ) = r] Xe ' Ps T' a+l3 . (133) 
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On the other hand, 



M{T a )M{T p ) = n x ^T' a+p . (134) 
Thus M preserves multiplication if x g -ph = x e -pf. This is indeed the case: 



i i \ j / \ k / 

where the last step follows from the fact that v T is the inverse of 7. 

We also need to show that for every a, M(T£) = [M(T , a )]L As in Ap- 
pendix B, we have 

M(Tl) = rf^tT'_ a (136) 

and 

[M(T a )]t = v q3 ' Ph TL a . (137) 

But the exponents are again equal because of the relation between 7 and v. 
Thus M must be conjugation by a unitary operator, as we wanted to show. 

2. Changing the value of w that relates the bases E and F 
The above result shows that a certain sort of change of basis preserves the 
translation operators up to a unitary transformation. In particular, this 
extends the result of part 3 of Appendix B at least to some other pairs of 
bases, but not to all allowed pairs. We also have to consider the case in 
which E is left unchanged and F is changed from (we~i, . . . ,we n ) to F' = 
(w'ei, . . . , w'e n ) where w' 7^ w. As we argue in Section 7, this change does 
not correspond to a unitary transformation of the translation operators, even 
up to a phase factor. Nevertheless, we can extend the result of part 3 of 
Appendix B to bases obtained by such changes. 

Let the translation operators T a be defined relative to the original bases 
E and F, and suppose it is known that for any unit-determinant linear trans- 
formation L, there exists a unitary Ul such that for any a 

U L T a Ul w T La . (138) 

Let T' a be the new translation operators, defined relative to the bases E and 
F' . From the definition of the translation operators it follows that 

T' a = T Ka , where K = f J Jj^ ) . (139) 
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We want to show that the new translation operators also have the above 
transformation property. 

To prove this, let L be any unit-determinant linear transformation, with 
Ul satisfying Eq. ([138)1 . Define a new unitary transformation Vl by the 
relation 

V L = U KLK - h (140) 
which exists since KLK~ X also has unit determinant. Then 

V L T' a Vl = Uklk^TkJJ^^ « T (KLK -i )Ka = T KLa = T' La . (141) 

Thus Vl is a unitary transformation that correctly transforms the translation 
operators, when the latter are defined relative to the new bases. This finally 
extends the result of part 3 of Appendix B to all pairs of bases that are 
allowed in our construction. 



Appendix D: Effect of linear transformations 
and translations on the indices a, 6, c, d, e 

1. Linear transformations 

Let L be a unit-determinant linear transformation, and let Ul be a corre- 
sponding unitary operator as in Eq. ([95)1 . Then, as we pointed out in Section 
6, if Q(X) is a quantum net, then 

Q'(A) = UlQ(LX)U L (142) 

is also a quantum net, covariant with respect to the same translation oper- 
ators. Moreover, starting with a given quantum net, one can generate its 
entire similarity class via the transformation (|142j) . with L ranging over the 
group of unit-determinant linear transformations. 

For the case N = 4, with our standard field bases E = F = (u, 1), we can 
label a quantum net by the indices a, b, c, d, and e of Fig. 5, which specify 
the quantum state assigned to each ray of phase space. Now, if we perform 
the operation of Eq. ()142|) on a quantum net, it is helpful to know how these 
indices change. In this section we present the relevant transformations of the 
indices for each of the three generators L 1; L 2 , and L 3 given in Appendix B. 
To define L3, we need to specify the field element z that appears in Eq. ([109)1 . 
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Let us choose z = uo. Starting with Eq. f|l 15|) . one obtains in the current 
setting the following unitary matrix to associate with L 3 : 

( 1 \ 

10 

1 

\ 1 J 



u. 



(143) 



The unitary operators associated with L\ and L2 are exactly as given in 
Eqs. (I113|) and (|114|) but specialized to the case of just two qubits. 

Let (a, b, c, d, e) be the indices characterizing the quantum net Q, and 
for a specific choice of L, let (a', b', d, d', d) characterize the quantum net Q' 
defined in Eq. ()142j) . From this equation and the matrices U{ given above, 
one can work out how the primed indices are related to the unprimed ones. 
Here we simply present the results. 
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2. Translations 

In the same spirit, we can consider translations of the 4x4 phase space and 
ask how they affect the indices that specify a quantum net. Given a quan- 
tum net Q(X) and a translation T a , consider the (equivalent but different) 
quantum net defined by 

Q'(A) = Q(T a X). (147) 
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Q and Q' can be specified by indices (a, b, c, d, e) and (a', b', c', d', e') as above, 
and one can ask how the two sets are related. Let us consider two basic 
translations, T(i,o) and T"(o,i)- One finds that for T(i,q) the primed indices are 
obtained from the unprimed ones by adding the vector (1, 0, 1, u, u>). In the 
case of 7( 0;1 ) the added vector is (0, 1, 1, 1, 1). (The latter result reflects the 
scheme by which we arranged the vectors in Fig. 5.) From these two cases 
one can obtain corresponding transformations for an arbitrary translation 
via a linear combination. 

3. Searching for invariants 

In Section 6 we introduced the function D(a,b,c,d,e) that identifies the 
similarity class of any quantum net for N — 4. To say that D has a constant 
value within each similarity class is the same as saying that it does not 
change when the quantum net is modified by either a translation or a unit- 
determinant linear transformation. Thus, one way to obtain the function 
D is to look for an invariant under all of the transformations given in the 
two preceding parts of this Appendix. One can show, in fact, that up to a 
constant factor and an additive term, D is the only second-degree polynomial 
in a, b, c, d, and e that is invariant in this sense. 
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